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We study the chiral fermion meson model which is the well known linear sigma model of Gell- 
Mann-and-Levy at finite temperature. A modified self-consistent resummation (MSCR) which re- 
sums higher order terms in the perturbative expansion is proposed. It is shown that with the MSCR 
the problem of tachyonic masses is solved, the renormalization of the gap equations is carried out 
and the Goldstone's theorem is verified. We also apply the method to investigate another known 
■ case at high temperature and compare with results found in the literature. 
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I. INTRODUCTION 



Several models have been proposed to describe hadron properties in the regime of low energies. Among these 
models, we adopt the linear a model of Gell-Mann and Levy Q which is a phenomenological model of Quantum 
Chromodynamics-QCD that incorporates two important features of QCD: chiral symmetry and partial conservation 
of axial vector current. The model was originally proposed as a model for strong interactions |l| , but nowadays it serves 
as an effective model for the low energy (low temperature) phase of QCD. It has the advantage of being renormalizable 
at zero Q and finite temperature j^] . Altough the linear sigma model lagrangian exhibits chiral symmetry, quantum 
effects break this symmetry spontaneously. Both from theoretical Q and experimental J5j point of view, there exist 
\ a great amount of interest in the study of chiral symmetry restoration at finite temperature. 

However, quantum field theory at high temperature has a well known problem that is the breakdown of the 
| perturbative expansion This happens in theories with spontaneous symmetry breaking (SSB) or in massless 

field theories because powers of the temperature can compensate for powers of the coupling constant. Resummation 
techniques which try consistently to take into account higher-loops are required. 

A systematic self-consistent approximation approaches based on the meson sector of the linear a model was previ- 
ously studied by Baym and Grinstein Jj|. After, Banerjee and Mallik |l(| proposed a modified perturbation expansion 
with the objective of calculating the two-point functions up to second order in the \<j) A theory. A resummed per- 
turbative expansion was proposed by Parwani fll~l| in order to go beyond leading order in the same model. More 
recently, Chiku and Hatsuda [Q in the study of the 0{N) 4 model presented a novel resummation adding a mass 
parameter determined later by the fastest apparent convergence (FAC) condition. We employ imaginary-time formu- 
lation (ITF) whereas in [ H real-time formulation (RTF) is used in the development of the optimized perturbation 
theory (OPT). In this paper we develop a modified self-consistent resummation at finite temperature and apply it to 
the investigation of the chiral fermion meson model. We study the temperature dependence of the chiral condensate 
and the effective meson and fermion masses by this self-consistent non-perturbative approximation up to one-loop 
order in the perturbative expansion. In the application of the MSCR to the study of the chiral fermion meson model 
at finite temperature, we divided the problem into three physical regions: low, intermediate and high temperatures. 
This is essential to identify the regions where resummation is crucial. In each region renormalization and satisfaction 
of Goldstone's theorem are discussed in detail. Our study addresses problems found in the context of the well studied 
0(4) linear sigma model and deals with a usually avoided point: the inclusion of the fermions. Also, we re-examine the 
chiral phase transition in static equilibrium in terms of the linear sigma model with our MSCR. Instead of demanding 
a infinite gap-equation, as has been done often in the recent literature, we perform the renormalization in stages in 
order to get finite gap-equations. 
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We also treat an explicit chiral symmetry breaking term in the Lagrangian which generates the realistic finite pion 
mass. Symmetry is never restored in this case. It is shown that in the limit of vanishing pion mass, namely when 
the chiral symmetry is exact, the inclusion of fermions does not change the order (nature) of the phase transition but 
only lowers the value of the critical temperature. 

This paper is organized as follows. In Section II we discuss the chiral fermion meson model and some of its features 
at zero temperature. In Section III the temperature is introduced via the partition function of the model which lead 
to the thermodynamical potential. The inclusion of loop corrections and the thermal gap equations are addressed to 
Section IV. In section V we apply the MSCR to the study of the massless A0 4 model in the weak coupling limit at 
high temperature. The renormalization of the self-energy is studied in Section VI. The numerical results are presented 
in Section VII. Section VIII is devoted to conclusions. 



II. THE CHIRAL FERMION MESON MODEL AT ZERO TEMPERATURE 



The Lagrangian density of the chiral fermion meson model which provides an explicit realization of chiral symmetry 
is given by 

C sym = ^ [i-fdp - g(*' + z 7 5 v? • r)] V + \ [(da') 2 + (dw) 2 ] - -^a' 2 + vr 2 - f 2 ) 2 , (1) 

where ip, a' , and tt represent the quark, sigma and pion fields, respectively, A and g are positive coupling constants 
and /,r is the pion decay constant in vacuum. 

If the up and down quark masses were zero, QCD would have a chiral SU(2)l x SU(2)r symmetry. In the vacuum 
this symmetry is spontaneously broken by quantum effects, with the result that there exists a triplet of Goldstone 
bosons. In reality the quark masses are very small but nonzero, so that chiral symmetry is only approximate and the 
pion has a small mass |l3| . An explicit chiral symmetry breaking term is added to the Lagrangian which generates 
the realistic finite pion mass so that 

with 

t-'symb — , (3) 

where c is small and positive. 

The term C sym is symmetric and invariant under an SU(2)l x SU(2)r chiral group and £ sy mb is the symmetry 
breaking term. Two Noether currents associated with (|l|), namely the vector current and the axial vector current, are 
given by 

-* — T 

= Ipjfi-lp + IT X d^TT, 

respectively. The equations of motion for the fields derived from the Lagrangian density ([!]) give the PCAC relations 

= or. (4) 

The effect of the term ctr' on the classical fundamental state, can be found by looking at the minimum of the 
potential 

V (*>,n) = ±(*> 2 +n 2 -f 2 ) 2 -c*i (5) 
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A(a' 2 +7? 2 -/ 2 K = 
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whose (unique) solutions are 



To first order in c, we have 



7? = 0, (8) 

ak 2 - / 2 k = c 



From (^), we see that a' has a non-zero vacuum expectation value. It is convenient to redefine the sigma field as 
a' — > a + v such that a has zero expectation value. As an effect of this shift the fermion field acquires a mass given 
by 

m.0 = gu (10) 
The shifted Lagrangian, C Sl of the new quantum theory reads 

C = -\(fl V 2 ? +*>- H" 3 v£ - - x )a + 
#7"^ - + ^[(dn) 2 - mln 2 + (da) 2 - m 2 cr 2 ] + 

-g~i>[o + i7 5 7? • t))iI) - -^[(t? 2 + a 2 ) 2 + 4i/cr(7? 2 + a 2 )] = 

-U(v) + C + £i (11) 

where U(y) is the mean field energy density, Co is the free Lagrangian and Ci is the interaction Lagrangian, defined 
by 

U{y) = ^{fl-u 2 f-cy, (12) 
Co = ^[i-fd M - + ^[(dn) 2 - mln 2 + (da) 2 - m\a\ (13) 

Ci = -g$[a + n 5 Tr ■ t)]i/> - ^[(tt 2 + a 2 ) 2 + 4va(n 2 + a 2 )}, (14) 

respectively. 

The meson masses read out of the shifted Lagrangian (flj) are 

ml=m 2 + \v 2 1 (15) 



m z a = m z + i\v z (16) 

where m 2 = —A/ 2 < 0. 

It is easy to see that the coefficient of the linear term in the sigma field, A(i^ 3 — vf 2 — j), in Lagrangian ( ) 
is identically zero by the minimal condition (||). This is due to the fact that the vacuum expectation value of the 
sigma field, (cr), should vanish at any order of perturbation theory fyj| , even if we include thermal corrections [^5| . 
The one-loop thermal tadpole corrections will modify this relation which will become temperature dependent. If v is 
allowed to be temperature dependent, the masses are temperature dependent as well. At any temperature, v is such 
that (cr) = 0. At zero temperature, when c continuously approaches zero, we have the solutions (7?) = and (a 1 ) = f„ 
which minimize the potential satisfying the Goldstone's theorem. 

The contact with phenomenology is made by fixing the parameters of the model to agree with the observable value 
of the particle masses in vacuum. Then, the tree level parameters of the Lagrangian are 

_ m 2 . - 3m 2 . 

A VI ' (17) 
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m v , ;0 (m^ ;0 -3m^ ;0 ) 
AK ;0 -2m^ ) ' liJj 

where m» ;0 = 139JWW , m CT;0 = 600MeV , my, ;0 = 340MeV and / w = 93Af e y. 

As we have mentioned earlier, our goal in this work is to study the chiral phase transition in the chiral fcrmion meson 
model and to analyze the thermal behavior of the temperature dependent meson condensate v and the meson and 
fermion masses. So it will be necessary to compute all the one-loop self-energies for the particles present in the model. 
Such self-energy diagrams have divergent pieces which must be renormalized if we want reliable results. In most of the 
approximations found in the literature [p|Jl^-[l8]| several difficulties have been found in the tentative of renormalizing 
the divergent gap-equations. Sometimes the undesirable parts have been ignored [fl9|| . The renormalization of the 
self-energy is studied in section|v| whereas the effective potential renormalization is performed in appendix ^[ For 
the purpose of renormalization it is necessary to add to C s a counterterm Lagrangian, Lc.t., needed to render the 
theory finite [§0), 

C = C s + Lc.t. (20) 

where 

Lc.t. = C.T. + D x m% + D 2 m* + D 3 mt (21) 

In (Ell) C.T. contains the appropriate counterterms to be used in the renormalization of the masses while 
-Di, 2,3771^, n a are necessary to keep the thermodynamical potential finite, as we will see further. As we are interested 
only in the study of the thermal effective masses at one loop order in the perturbative expansion, other counterterms 
necessary to renormalize the coupling constants are not explicitly shown. 

III. THE PARTITION FUNCTION OF THE MODEL AND THE LINK TO STATISTICAL MECHANICS 

One of the most fundamental objects in thermodynamics is the partition function, defined by, 

Z = Tr e-P*, (22) 

where H is the Hamiltonian of the system, (3 = 1/fcgT = T _1 with the Boltzmann constant, fc^, set equal to one and 
the trace, Tr, in eq. ( ^2| ) meaning the sum of the elements of the matrix e~^ H in all independent states of the system. 
All information concerning the equilibrium thermodynamic macroscopic properties of the system are obtained from 
Z. 

In relativistic quantum field theory, the partition function can be derived from the Feynman's functional formalism 
plf . The bridge between quantum mechanics and statistical mechanics is achieved by the heuristic introduction of a 
variable defined as r = it. Also, the fields are constrained to obey periodic(anti) boundary conditions: <fi(r, 0) = <f»(r, 0) 
for bosons and •(/'(r, 0) = —ij){r,[3) for fermions. Following these prescriptions, we get 
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Z[%j), ip, cr, tt] = N' I V[(j>] exp[ / dr \ d 3 xC{<j), d<j))] = (23) 
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Here we have introduced a short hand notation for the Euclidean space-time integral: S — f„ d 4 xC = dr J d 3 xC, 
T>[4>] is an abbreviation for the integral over -0, ip, a and 7r, N' is an unimportant infinit constant and C is given by 
(§)• 

Next we introduce the thermodynamical potential, f2, dinned by 

n(2» = ~ln3 (24) 
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where InZ — ^(U(u) + Di^^m^ n a ) + IuZq + InZj. Since the interaction action Sj contains terms which are more 
than quadratic in the fields, it is not possible to carry out the functional integration above in closed form. For a while 
we will neglect Si in our calculations. This amounts to considering only the tadpole contributions. Thus, 



d 3 p 
(2tt) 3 



— u n 
2 



Tln(l 



fii(T, v) = \{fl~ v 2 f - «/ - [D x m% + D 2 mi + D 3 m 4 a ] - ^lnZ = 

- v 2 f -cv- [D im % + D 2 mi + D 3 m 4 a ] + 



(25) 



+ 3Tln(l - e - ^) - 2 • 2 [uty + 2Tln(l 



with uj 2 a = p 2 + to 2 , w 2 = p 2 + to 2 and w 2 = p 2 + m 2 . In the third line of ( p5| ) the first factor 2 multiplying the 
bracket which contains the fermion contribution, comes from the spin degrees of freedom, whereas the other factor 
two is due the isospin degrees of freedom. Inside this same bracket there is another factor 2 corresponding to the 
particle and antiparticle contributions. The thermodynamical potential f^i is precisely the one-loop effective potential 
|^|,^) of the linear sigma model, and it can be expressed as, 

Q^T, v) = U{y) + D^m 4 ,^ + Q?(z/) + Slf{T, v) (26) 

The equation of state of the noninteracting system composed by a (free) relativistic boson and fermion gas is 



-T 



d 3 p 



P = ^\nZ = -^(T,p) = 
[ln(l - e~ pu ') + 31n(l - e'^) - 81n(l + e~ 0u *)] 



(27) 



where Pq is the thermal pressure. 

The integration over the temperature independent terms f2?(j/) = J [f^g + f^w ~ 4w^] actually diverges. It is 



exactly the counterterms Dx^rrA v a which will take care of these divergences. We discuss this further in section VI 



It can be known, from thermodynamical considerations, that in thermal equilibrium Q is a minimum with respect to 
variations in v p3| . Applying this extremum condition, we have 



gOi(T, v) 
dv 



= 



gQiM _ r d 3 P 



dv 



/ts& 3A(^ + ^)-8 3 



2 1 



(28) 

j/ represents the sum of 
l 



Since Q^(T = 0,v) = 0, the divergent quantity 

the tadpoles at T = 0. On the other hand, ^hSE^l [ s finite at T ^ because of the natural cutoff 

present in the integrals. We could let the counterterms absorb the finite parts of dC g^ together with the infinities 
since the vacuum contribution is irrelevant for the discussions of thermodynamics. But, as we will see in section IV C , 
these terms are important in the verification of the Goldstone theorem in our self-consistent treatment. After these 
considerations, equation (|2q) is written as 



where the functions F and G arc defined as 
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u^ieP"* + 1) 



(29) 



(30) 



3A f TYt^ 771" \ TH 

G {un^im*,^)) = j—^ [ml M~^t) + m l ln ("7/) ) ~ J^2 m l ln (~^)' ( 31 ) 

respectively. In eq.(|3l|), /x is the renormalization scale. 

As a first approximation, we consider only the thermal loop corrections to the effective potential. This approximation 
allows us to get an analytic expression for the approximate critical temperature. 
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d 3 p 

(2^)3 



v - X(f 2 - v 2 )v - c = 



U>x U! a J City 

where n,^ and are the usual distribution functions for bosons and fermions given by 



1 



(32) 



(33) 
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(34) 



respectively. In the above expression, when one minimizes the effective potential, one is summing the thermal tadpole 
contributions to the usually called mean field equation. The chiral condensate, v 1 which is a non trivial solution of 
this integral equation now depends on T. This equation can be solved with an explicit analytic form in the high 
temperature limit. The leading terms in the high temperature approximation for this integral equation are 
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(35) 



The above equation has a real solution that is a slowly decreasing function of temperature, but does not vanish. Thus, 
when c ^ the symmetry is never restored. On the other hand, when c = the non trivial solution of (B5T) is 



* 2 = />-i(i 



— )T 2 
3A ' 



The critical temperature is defined as the temperature where the condensate goes to zero. It is given by 



T . 



(36) 



(37) 



It shows that the inclusion of fermions does not change the order of the phase transition, but only lowers the value 
of T c . Note that the interactions of the mesons with the fermions forces the "critical" temperature to depend on the 
coupling constants A and g. If g = we recover the result of p2|. 



IV. INCLUSION OF LOOP CORRECTIONS 
A. The first necessity: Beyond the mean field approximation 

Let us analyze the finite temperature behavior of the tree-level meson masses ( |l5|) and (|l6| ) as functions of the 
thermal expectation value of the sigma field, v(T). Since v 2 {T) decreases as T increases and m 2 < 0, the particle 
masses becomes tachyonic. Another problem which arises is the fact that Goldstone's theorem is not satisfied in the 
ordered phase (when c = 0), i.e., replacing ( |36"| ) on (|l5|), we obtain a non zero pion mass given by m 2 = — -|(l-|-^-)T 2 . 
This pathological behavior is due to the fact that in our approximation we have neglect the interaction action Si 
in the thermodynamical potential (^4|). The result is that the mean field approximation can be trusted only in the 
approximate prediction of a phase transition at T c given by (|37]) . It is incorrect in what concern the description of the 
finite temperature behavior of the meson and fermion masses. So it is necessary to include all one-loop corrections 
from all 1PI diagrams present in 57 to the masses. 

Following the program of [fl3[ we expand the partition function in powers of the interaction, in order to get the 
one-loop self-energy corrections. 

lnZ,-ln(l + ^- ) (38) 

n— 1 L J 

The one-loop 1PI graphs come from InZi + lnZ2, which are given by 



J [d(/)}e So Si 



1^1 = (39) 



G 



1 JMe*>S l 2 l JMeSoSj 
lnZ *--2 { \dd>]eSo } + 2 \d6] e s° (4 ° } 



where the disconnected diagrams cancel in InZ^ and the diagrams which gives rise to tadpoles in the self-energy are 
not to be included, since their effect is already considered in the mean field equation. The terms which "survive" 
come from: 

lnZ 2- l0 o P = X j[d^oJdrJd 3 X (^+a^ 
1 4 f{d(j)}e s ° 

2 X dTldT2 J d Xld X2 JW^ 

The 1PI graphs from this expression can be represented diagrammatically as shown in Fig.|l|. 
The self-energy for bosons and fermions are defined, respectively by 

D^k)"^ =D a, OT K ) k)- 1 +n (T , 7r (w n ,k) (42) 

©K.k)" 1 =D K,k)- 1 +EK,k) (43) 
where Doo- !7r (w„, p) and 2?o(wn,p) are the tree- level boson and fermion propagators, expressed respectively as 

Do^Wn.k) -1 = cj 2 + k 2 + m 2 jW (44) 

PoK.kj-^^F-m^ (45) 

Here, to n are the Matsubara frequencies, defined as uj n — 2nnT for bosons and uj n = (2n + 1)ttT = u> n f for fermions. 
To one-loop order the self-energy expressions are given [fl3f by 

Ci y2 — loop 

IW = -2(^ W (46) 

^(^-f^, (47) 

The self-energy graphs to each particle can be pictorially represented as by cutting one of the corresponding loops 
in the diagrams representing (In Zj~ oop )ipi. After the integration in x, r and in the fields in (In Z\~ lo ° v )\pi and the 
differentiations above, we obtain the following expressions for the self-energies at one loop order 

lMfc ,k) = ]Tn„ = A5T^ [ ^D .K,p) + AT^ f ^D 0CT K,p)+ (48) 

n a (fc ,k) = ^n CTs: = 3AT^T / ^_D 0ff (o; n ,p) + 3AT^ / ^D^^.p) + (49) 

-6\ 2 is 2 Tj2 J ^3D 07r (w„+i,p + k)D 07r K,p)-18A 2 ^ 2 r^ J ^^D 0<7 (w B +,,p + k)D 0ff (w n ,p) + 

d 3 p 



T Y, J 7^' & [ Z? 0^(w B +I,P + k)I>0*(Wn,P)] 
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E(fc ,k) = ^s l = -. 9 2 r^ / 



d 3 p 



(27T) 



3 D 0ct {uJn+l , P + k)2? ^ (W„ , p) + 



(50) 



j—^ D 07r (w n+; , p + k)D </> (w„ , p) 



The diagrams representing the pion, sigma and nucleon one-loop self-energies are drawn in figures |,| and |. 

We note here that we could get the same results for the self-energies directly applying Feynman rules to construct 
the diagrams with the appropriate substitutions: the (^-function at each vertex is replaced with a Kronecker delta 
which imposes conservation of the discrete energy (fcg = iw n ), and round each loop of a thermal graph with J 



which are the finite-temperature Feynman rules p3 23 1. When the summation over n is performed, each 



(3 I (2tt) 3 

graph of the self-energies is separated into two parts, namely a temperature independent part (at this stage), which is 
divergent, and a temperature dependent part containing the Bose-Einstein distribution in the case of bosons or the 
Fermi-Dirac distribution for fermions. 

We will adopt the definition of mass at finite temperature as the real part of the pole of the corrected propagator 
at zero momentum (k = 0). Thus, from eqs. (ff2|) and (43) we have 



d^k, |k| = or 1 = cui 



-k 2 



• n»( Wni |k| = o) = o -> 
•n*(T, VrJki = o) = o 



(51) 



o)- 



-k 2 



- n CT (w n , |k| 
n CT (T, /cq j(T . 



= 0) = 
|k| = 0) = 



(52) 



P*W, |k| = 0)- 1 = 7^ - m v , + S a 



7 M S M = 



(53) 



fe ,v> + ^o(T, korf, |k| = 0) + S S (T, ko,,/, |k| = 0) - m$ = 0, 



where the arrow indicates an analytical continuation from discrete to continuous energies in Minkowski space. Hence 
the physical masses are the values of the k^^^ which are the zeros of the functions (|5l|), (52) and ( |53| ) above, i.e., 
the location of the poles in the limit k = 0. The full self-energy expressions 11^, IT-, T,q and E s are shown explicitly 
in appendix [B[ The renormalization of the self-energy is studied in section VI . Through out this paper, we will 



use dimensional regularization, but omitting, for notational simplicity the factor /i 2e which multiplies A. Since our 
calculations does not require traces involving an odd number of 7 5 matrices, we use the definition of j 5 as in p4| , p5"| . 

Since these expressions are self-consistent they have to be solved numerically. For each fixed temperature one finds 
a value of M which satisfies the equations above. On the other hand, if one is interested only in the meson sector 
of the linear sigma model, the integrals in the self-energies could be evaluated exactly in the high temperature limit 
and at low frequency where the boson diagrams involving three-point vertices which are proportional to \ 2 v 2 may be 
neglected. This is not consistent if one wants to study the behavior of the condensate v and the particle masses in all 
ranges of temperatures. It is important to note that when c = the three-point vertex boson diagrams are significant 
in the region Tj(~ 0) < T < 7>(~ v(T)) and when c ^ 0, v{T) ^ for any finite value of T. 



B. The second necessity: The resummation 



The expressions for the self-energies appearing in eq. (pl[), (|5^ ) and (53) are functions of w I)(7 ^ which are expressed 
in terms of the mean-field masses. As we discussed, the meson masses become negative as the temperature increases. 
Thus, in the computation of the one-loop corrections, the masses running in the loops become tachyonic. A proper 
resummation of higher order loops is naturally necessary H . Various resummation methods have been proposed in a 
tentative of curing the problem of the breaking down of the perturbative expansion at high temperature. In effective 
models, when a phase transition occurs, one can find tachyonic masses even below T c . The O(N) linear a model 
which is one of the laboratory effective models employed to study QCD has been investigated by different authors 
using different techniques. One of these methods is the CJT formalism p(| which provides for a consistent loop 
expansion of the effective potential in terms of the full propagator. The CJT formalism elegantly provides for a gap 
equation from stationarity conditions for the daisy and super-daisy effective potential. However some authors use 
this non-perturbative approach with an ansatz for the full (corrected) propagator in which the thermal corrections 
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are momentum independent. These corrections are the finite piece of the divergent integral (which is temperature 
dependent through the gap equation for M) plus the finite explicit temperature dependent piece. This is the Hartree 
approximation, and means resuming only the "bubble diagrams" that are dominant at high temperatures. Another 
non-perturbative approach widely found in the literature is the large- N approximation. The N — > oo limit facilitates 
the calculations, but can lead to inaccuracies |^,^9). One must be careful in taking the large- N limit since its 
truncation depends on the problem to be studied and the relevant value of N [^7| . In this case the three-point vertex 
diagrams are omitted which in principle makes sense only in the N — > oo limit since these sunset diagrams are of order 
1/N. In these two kinds of treatment one can not study the bosons interacting with fermions (with the interactions of 
the linear sigma model) since the self-energy diagrams are momentum dependent which invalidates the ansatz cited 
above. It is worth to remember that the Hartree approximation does not satisfies the Goldstone theorem p9| , p8[ . 
This fact may be attributed to the non inclusion of these diagrams. Although the finite temperature mass in these 
approaches is the pole of the corrected propagator, it is not the true mass, since their Il„ i<T are not the true one-loop 
self-energy functions (see discussion below). The corrections included only shift the masses. Once we are interested 
in the study of the masses behavior also in the range < T < T c (if c = 0) the three-particle vertex diagrams will not 
be neglected. The inclusion of these diagrams brings an additional complication since the self-energy now depends on 
the momentum. 



C. A non-perturbative resummation method: The MSCR 

Let us now introduce our procedure which resumms higher loop diagrams in the mean-field (tree-level) propagators. 
The method consist in recalculating the self-energy, in steps, using in each step the masses obtained in the previous 
one such that M 2 = (A n + 1)M^_ 1 + II(M n _i), where n is the order of the non-perturbative correction and A n is 
the coefficient of the appropriate counterterm. With this procedure it is easier to identify and absorb the divergent 
parts of the self-energy in order to have finite gap-equations. The goal is to make renormalization possible since the 
masses which multiply the divergences are necessarily the same as in counterterms. 

Application of MSCR 

Here we apply the MSCR in the study of the chiral fermion meson model at finite temperature. The analysis of 
the problem has to be done carefully which will be divided into three regions. 

Region I: The low temperature region 

The first region is for < T < T*, where T* is the temperature where v(T*) — /„.. This implies that m 2 = and 
consequently the appearance of infrared divergences in the self-energy. So, in this region 
Step 1: 

Start with the mean-field effective Lagrangian where the condensate and the masses are given by: 

[F (^^(m^^o)) T) + G (^^(m^^o)) - A(/ 2 - v 2 )} v - c = 0, (54) 

M?, = ml=m 2 + Xv 2 , (55) 

M 2 = m 2 a = m 2 + 3 At/ 2 , (56) 

M^,fi = m^j = gv. (57) 

Step 2: 

Evaluate the one- loop self-energy corrections to these masses from the equations presented in appendix [b| and define 
the condensate and the first order corrected masses as 

[F K^M^.o^T) + G K^M^^o)) - A(/ 2 - v 2 )} v-c = 0, (58) 

M 2 A = Ml, + UAMl Ql M 2 , M^o) = (59) 
{A l +E X + l)m 2 + {A x +E X + l)\u 2 + II*(M£ 0) M 2 , M^ ) = 

m 2 + A^ 2 + nf "(Ml o, M 2 , M^o), 
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Ki 



{B x + Fx + l)m 2 + (Fx + l)\is 2 + (Ex 



= Ko+n*(Ka>Ka> M *fi) 
l)2A^ 2 + n CT (M 2 . ,M 2 ,M^ ) 



3A^ 2 



(60) 



n 



Ren 



M^x = M^ Q - E(Af£ , Af£ , Afy,o) = 
(Cx + l)gu - Z(Ml , Ml , Afy,„) = 
^ Ren (M 2 fi ,M 2 fi ,M^ ). 



(61) 



where IL, 



M v a \ k 



0)° + ^(T; ko. 



0;k = O)' 3 , Ai, Ax, B t , Bx, Cx, Ex and Fx are the 



appropriate coefficients of the counterterms added to the mean-field effective Lagrangian needed to render the model 
finite up to this order, which are shown in section VI. The requirement that fco l7r .<T = in 14^ a excludes the possibility 
of tachyonic tree-level masses since thermal effects provides for the pion a non-zero width due to the Landau damping 
process. For the fermions we adopt the requirement that £ — S(fco lV > = M^;k = 0)° + S(T;fco j ^ = 0;k = 0)^ to 
prevent a similar consequence. The resummation has to be done exactly to avoid this problem. In this range there is 
no necessity of resummation since the masses running in the loops are positive, i.e., M 2 > 0, M 2 > and > 0. 
Renormalization: 

The renormalization is done normally since the masses multiplying the divergences are the same as in the counterterms. 
The coefficients of the counterterms are found in section VI. 
Goldstone's Theorem: 

In the exact chiral limit (c = 0) and low temperature phase (where v ^ 0) from eq. ( |5l| ) at (ko — > 0, |k| = 0), we 
have 



Mix = Ko 



Koto 



\ 2 v 2 

' M 2 , ~ M 2 , 



-4 

(4tt) 2 



(4tt 

Ml^{%)-Ml ln(^ 





-m 2 + \iy 2 + U Ren (k -»0,|k| = 0) = 
- (4^)2 M <W>H 



(62) 



5A m2 . (M 2 A 5A 
^y^oln( — )+y 



W<r(M ff>0 ) 

(47r) z ' \ efi 




g 2 ^(%) 

cj^(M0 >o ). 



- 2AV 



dpp 2 
~^ 2_ 




(47r)^ ^ \ ep^ 



dpp 2 rty = 

7T 2 City 



since M 2 Q — M 2 = 2\v 2 (in deriving eq. ( |62| ) we have used eq. ( p9 



Region II: The intermediate temperature region 

This is the region for T* < T < T c where the resummation is really necessary, since the masses in the loops 
are zero or tachyonic. The problem here is more complicated than in the other temperature regions since phase 
transition takes place. As is well known, around the critical temperature quantum fluctuations become essential and 
one loop corrections may not be enough. In fact, as we proceed to show here, this scheme (i.e., allowing only one 
loop corrections in the perturbative expansion) prevents us from adequately renormalizing the gap equations for the 
tree-level resummed masses and satisfying Goldstone's theorem. 

We show next that this is indeed the case and that a possible (but inconsistent !) way of achieving renormalization, 
would be e.g. to consider only the diagrams for which IT^ = H a . This is very frequently implemented in the literature 
H,HJl9|j2f . 

We will adopt the point of view that the MSCR fails in this region and higher order loop corrections in the 
perturbative expansion will be necessary at this stage. This is subject of current investigation. However we explicitly 
show were the problem appears. 
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Step 1: 

Start with the mean-field effective Lagrangian where the condensate and the masses are given by: 

[F (ci} Vt<r ,^(m %!l j^-o),T) + G (w 7r ,a,v>(nv,c7,V'io)) - K.fl - v2 )\ u-c=0, (63) 

M* „ = ml = m 2 + Xu 2 , (64) 

M 2 = ml = m 2 + 3 \v 2 , (65) 

M^,fi = m^j = gv. (66) 

Step 2: 

Evaluate the one- loop self-energy corrections to these masses from the equations presented in appendix [b| and define 
the condensate and the first order corrected masses as 

[F (^^(M^^o),! 1 ) + G (wtt, 0-^(^,0-^,0)) - Kfl ~ v 2 )] v — c = 0, (67) 

Ml, = M 2 >0 + H„(M% )0 , M 2 , M^.o) = (68) 
(A, + E l + l)m 2 + (A, + E,+ 1)A^ 2 + n*(Af£ , M 2 , M^ Q ) = 

m 2 + \v 2 + nf "(M 2 j0 , M 2 , M^o), 

Ml, = Ml + n a (Ml , M 2 , M A0 ) = (69) 
(B, + Fi + l)m 2 + (Fx + 1)A^ 2 + (B, + \)2\v 2 + n CT (M 2 i0 , M 2 , M 0jO ) = 



111 



1 < 3\v 2 +n*™(M 2 Ml „,M^o), 



Af^,i = M^, - S(M 2 , , M 2 o, M^, ) = (70) 
(Ci + - S(M 2 i0 , M 2 , M^, ) = 
gv-V R ™(Ml ,Ml ,M^ ). 

At this stage of the application of the method, the renormalization is possible, Goldstone's theorem is verified, but 
the tree-level masses are zero or tachyonic. 
Step 3: 

Now we take the masses computed in the previous step and improve the results defining a next-order non- 
perturbative correction. With this we get a new effective Lagrangian where the condensate and masses are given 
by 

[F T) + g {cj^AM^a)) - Kfl - ^ 2 )] v-c = 0, (71) 

Kfl = Ml, + IL^Ml,, Ml „ M^) = (72) 
(A 2 + E 2 + l)m 2 + (A 2 + E 2 + l)\v 2 + 
(% + E 2 + l)nf "(M 2 i0 , M 2 , M^o) + n T (M 2 4 , Ml „ M^,) = 

m 2 + \v 2 + Hf "(M 2 !l; M 2 17 M^) - T^nf "(M 2 il7 M 2 x , M^,)l, 

M 2 2 = Ml, + U a (Ml„Ml„M^ A ) = (73) 
(B 2 + F 2 + l)m 2 + (F 2 + 1)A^ 2 + (fl 2 + 1)2A^ 2 + 
(I 2 + F 2 + l)Hf "(M 2 i0 , M 2 , M^,o) + rLfM^M^Mw) = 

3A 1 
m 2 + 3A^ 2 + U Ren (M 2 ,,M 2 ,, M^,) - -—-Il Ren (M 2 ,,M 2 ^,,M 1 f,,,)-, 
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(C 2 + \)gv - (C 2 + 1)E 



/)'r/! 



A^,2 = M^,i 
(M» 0s M» 0s A^, o ) 



E(M» 1S M» 1) M^, 1 ) = 



(74) 



It is shown that in the A</> 4 model, this (first) recalculation is equivalent to the sum of an infinite set of diagrams, 
namely the "daisy" sum || or the set of ring ]l3|] . The coefficients of the temperature dependent mass countert- 
erms A 2 , B 2 and C 2 are fixed in a manner to cancel not only divergences proportional to H Ren (M 2 , M 2 , Myo), 
H Ren (M 2 Q, M 2 o, M^fl) and T, Ren (M 2 , M 2 , M^ jQ ) respectively, but also these terms together. That is, at each 
stage of the procedure, for n > 1, in the expressions for M^^^^n, the self-energy n ffi(r (M 7r)tri ^ in _2) (or S(M WiCTj ^, jrl _2)) 
have to be cancelled to avoid overcounting of diagrams. 

This shows explicitly that renormalization can not be performed within this approximation scheme. This is not 
surprising since in this temperature region quantum fluctuation may need a more thorough description. So, there is 
no reason to believe that only the "daisy" diagrams should be resummed at low and intermediate temperatures. In 
fact, the "daisy" graphs contributions are dominant at high temperature |q]. 

Renormalization: 

Since in this region M 2 n — M 2 n = 2A^ 2 + All, where All = IIo- — H^, there is the presence of undesirable non- 
rcnormaliziblc terms. These terms are the last ones on the r.h.s. of equations ( |72] ) and d73| ) which come from equations 
(B2) and ( |B(i| ) respectively and can not be absorbed in the counterterms. 
Goldstone's Theorem: 

In this region, Goldstone's theorem is satisfied only if All — > and the contribution in (B3) is decoupled into integrals 
proportional to A. This would assure the cancellation of M w at (fco — * 0, |k| = 0). The reason for this frustration is 
the same as for the lack of renormalizability. 

Step 4: (applicable only in the case where All = 0. This guarantees that renormalization and Goldstone theorem 
can be satisfactorily implemented at each step. In particular this will be the case for the high temperature region as 
will be shown next.) 

Proceeding with the iteration, in the limit n — > oo the masses M n have formally the same expressions as the masses 
M n _i which are already renormalized. Thus, in this limit we will have, 



[f + G (wi,^(Af,,^,„)) - A(/ 2 - v 2 )\ v - c = 0, 



(75) 



(76) 



3A;/ 



n ,M* My, n ), 



(77) 



M^ n = gu-E Ren (Ml n ,M* 



(78) 



At each intermediate step, in the loops we set k 2 n ^ = M 2 . n _ 1 in the computation of M 2 g , n . This ensures 
the cancellation of the divergences in all stages of the process, since the masses in the counterterms will necessarily 
be the same as in the divergences. In the end, in the resulting equations of interest (to be solved nummericaly) , 



K 2 



M 2 ^ as it should. By our MSCR we have gotten a set of four coupled non-linear integral equations to be 



solved self-consistently, with finite gap equations for the tree-level masses, which read 

Ml 



m 2 + \v 2 + IL Ren (fe 0w = M w , |k| = 0) 
M 2 =m 2 + 3\v 2 + U Ren (k 0n = M n , |k| = 0) 
My =gu- X Ren (k y = My, |k| = 0) - Sf n {k y = My, |k| = 0) 



(79) 
(80) 
(81) 



m 2 + A^ 2 + 



3A 
(4tt) : 



dpp 2 



3A 



n a {M a ) 
u a {M a ) 







8 5 2 






(4tt) 2 


n-xiMn) " 


) +8.g 







2 ny(My) 

LOy(My)_ 



(82) 
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where the renormalization scale fi can be determined by a physical condition. We choose fi such that the pion mass 
has the correct value at T = 0. 

Now we have to go back to real world encoded by eq. (|l| to |53|) . In order to get finite physical masses from the 
pole of these equations it is necessary to sum and subtract the finite quantities n^ en and £ flen which will be regarded 
as mass parameters J7|,^0|-p^] . This corresponds to the reorganization of the perturbative expansion. Now we rewrite 
eq. @ as 

C = -hfl - v 2 ) 2 + az + ^iyfy - + 

4: 

i[(<97?) 2 - M 2 n 2 + (da) 2 - M 2 cr 2 } - g~$[a + ij 5 if • r)]V> + 
~[(t? 2 + a 2 ) 2 + 4va(n 2 + a 2 )} + S flen ^V + "(t? 2 + o 2 ) + C.T. (83) 



The last two terms on the third line of cq.(jS3D must be considered as extra "interaction" terms and will naturally be 
present in \nZj, eq. ([4l|) . The counterterm structure of eq.(|83]) is the same as the one present in the method for the 
pion and fermion, differing only by numerical factors in the case of the sigma mass renormalization. By eqs.(^) and 
([|7]) the extra contribution to the self-energy read 



Yjextra jjextra j-|-.Ren (84) 



where we have defined 



n 



Ren 



n ff (fc , w = M^k = o)° + n ff (t, fc - o, k = o)' 



^Ren 



E(fc = M^,,k = 0)° + S(T, k = 0, k = 0) / 



(85) 

(86) 
(87) 



As a result, the final resummmed tree-level masses (eqs.(|79|), ( |80| ) and (|81| )), may be used in eqs. (|5l| ), ( p2| ) and (|5J 
if one wants to study for instance spectral functions as the authors of |12j ], or decay width as in |p8|. 
Goldstone's Theorem: 

If the sunset type graphs are neglected, the self-energy function at one-loop is not complete, despite of the fact that 
the definition of masses as poles of propagators at zero momentum is still valid. Now, having the result for the full 
one-loop (and higher-order loops contributions from the resumation) self-energy function we can test algebraically the 
fulfillment of Goldstone's theorem in the exact chiral limit (c = 0) and low temperature phase (where v 7^ 0). From 
0, |k| = 0), we have 



eq. @ at (fc 
Ml 



n 



total 



M 2 + u v (k ^o,\k\ = o) + n* 



J 3A 






\(4tt) 2 





'Ml 

ep 2 



Mi In 



m 2 + \v 2 + n 
M 2N1 



(k 



dpp 2 







3A 



5A 



A 2 *, 2 M < 



(4tt) 2 



M 2 
ep? 



Mlln ^ +^ 



5A 
~2 



-M 2 ln(^ 



~(4tt) 2 



Ml - Ml 



2tt 2 
50 dpp 2 n v 

- 2AV 



n a (M a ) rhrjMn) 
Lo<j(M a ) w,(Mj) 

A , (Ml 
V j-r^Ml In — I 

(47r) 2 a \ep 2 




W^(M^)_ 

A f 00 dpp 2 n 



dpp 



2 r 



1 



tin 



1 



Ml - Ml 



% 2 Af2l 



K 

ep? 



uj a Ml - Ml 

/2 f°° d PP 2 n i> 
Jo 



= 



7T^ 



since M 2 — M 2 = 2A^ 2 (in deriving eq.(p8[) we have used eq.(p2|)). 

To our opinion the fullfilment of Goldstone's theorem is ultimately related to the preservation of the relation imposed 
by chiral symmetry to the tree-level masses. Moreover, it is crucial to keep all diagrams of a given order. This is due 
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to the fact that, strictly speaking, a loop expansion is an expansion in powers of the Lagrangian. As discussed in |l5| ] 
in order to respect the symmetries of the Lagrangian, one must retain all diagrams to the given number of loops. 

Then, in the absence of the explicit chiral symmetry breaking term, one has, 
for < T < T c 



Ml 



Mi 





2\v 2 



(89) 



for T = T r 



Mi = Mi 



(90) 



and for T > T c 



Mi = Mi 



3A 



dpp 2 n b 

7T 2 LU b 



•V 



dpp 2 rif 

IT 2 Uf 



(91) 



which shows chiral symmetry restoration. Here b stands for bosons and f for fermions. We could interpret the result 
in the r.h.s. of eq. (|9l| ) as if each independent pion effectively "sees" one sigma and the other two pions and four 
fermions (since /i, the chemical potential, here is zero). On the other hand the sigma "sees" the three pions and four 
fermions. This equation serves to define the critical temperature in which the common masses of the particles vanish. 

9 2 f 2 

In the high-temperature limit of these integrals, we find that Tf — — 
eq.f 



2g T- , as predicted by the mean-field analysis in 



Region III: The high temperature region 

This is the region of high temperatures, T > T c , if c = and v = or T > Tj, where Ti is defined as a "inflexion" 

,2 



temperature, for the case c^O and v « f„ such that M 2 

n.(Af w 2 



A£ 



Mi , 

7T, 1 



2 0l M 2 M^ ) 



m 2 + n^{Ml fi ,Ml fi ,M^ fi ), 



(92) 



rn 



■ 2 +U a (Ml ,A4l 
1 Ml fi ,M ilfi )=M\ 



m 



,0) 



M^o) = 
M^ ) = 
\i = Ml 



(93) 



M^,i = M^.o - X(Ml ,Ml ,M^ i0 ) = 
gy-E(Ml , M* M^, ). 



(94) 



Note that the pion and sigma masses become degenerate and the problem encountered in the previous region 
(T* < T < T c ) is no longer here since ALT = in this region of temperatures. In this case, the masses in the loops 
can be neglected, and we have 



M 2 = (A 1 + l)M 2 + n(Afo) = m 2 + 



A 



1 



2 j1 

3A 



(95) 



11 
T 2 



If we set 3 = these results agrees with the ones obtained by Bochkarev and Kapusta (22| 
Following the iterations, we find for the n-th iterated mass 



Mi 



(A n + l)M^_ 1 +n(M rl _ 1 ) 



(96) 



3A 



1 - -M B _! 
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In the limit n — > oo, we get 



which can be easily solved for M, 



M z = mf + - I 1 



3A 



)t 2 


l- A m 







(97) 



M 



For T >> T c 



3Xf « T Y + Xf*( T2 1 



M2 = l(^ Tl 



3A/ 2 T 
2^ T 2 ' 



(98) 



(99) 



V. THE MASSLESS A<^ 4 AT HIGH-TEMPERATURE 



Now we apply the MSCR to study a very popular model: the massless X(f> 4 model in the weak coupling limit 



(100) 



Ml = M 2 + H(M ) = 



Mo = 
AT 2 
~24~ 



(101) 



at this stage of the procedure there is no necessity of adding countcrterms since up to this order there are no ultraviolet 
divergences in dimensional regularization fTjj] . Here II is the 1PI one- loop self-energy to lowest order, namely the 
"bubble" of Fig. |(a). 



Ml = M 2 + n(Mi) 



= (A 2 
AT 2 
24 



l)U Ren {M ) + n(Mi) 
3Mi 



(102) 



1 



ttT 

A 
24tt 2 



0(A 2 In A) = 
+ 0(A 2 In A) 



with the result that this correction to the mass is of order A^, which is an signature of the non-perturbative re- 
summation. The temperature dependent counterterm is fixed so as to cancel the divergence and avoid overcounting 
of diagrams, as explained before. So, A2 = — 1 + w^Fp e* ^ e diagrams used, in a given number of loops, in any 
resummation method must be the same in all stages of the process. What changes is the masses running in the loops at 
each iteration. This is because one must keep the same fundamental theory in the recalculation of the self-energy. The 
result shown in eq. (102) is in agreement with the one obtained by Parwani's resummed perturbative expansion |Tl] ] 
(see eq.(2.12) of his paper). The second iteration corrected mass, M2, which was obtained in our method evaluating 
Fig. ||(a) with Mi in that loop can equivalently be achieved calculating the "daisy" sum, that is a summation of 
the infinite set of "daisy" diagrams of Fig. |](b) with Mq in the loops. In this case all "daisy" types diagrams are 
IR-divergent since M Q — 0, but their sum is IR-fmite Jl^Jll[| . 
Continuing the iterations, we find for the next correction 



Ml = Ml + n(M 2 ) 



AT 2 
~2T 



'3 

3Mi 



(103) 



3Mi 



When A < < 1 we get 
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and for the n-th iteration, we obtain 



AT 



3=1 



2J- 



/ A 



V247T 2 



(105) 



The "superdaisy" sum || corresponds to the limit n — > oo of equation (105) and it can be summed up (for A << 1) 
to give 



M 



AT 2 
~24~ 



1 



3 f A \ 2 
2 V24ir 2 / 



1 + 3 ( A \2 
. 1 ^ 2 V24^V 



(106) 



VI. RENORMALIZATION 



A. Determination of the counter-terms 

The divergences are regulated via dimensional rugularization. To renormalize the divergences, we use the Minimal 
Subtraction scheme where only the poles are eliminated by the appropriate counterterms. The first-order parameters 
of the temperature-dependent counterterms read 



6A 1 -r 



{At:) 2 r 



12A 1 

:~, -&1 



(4tt) 2 T 



with i 



4-d 



7 + log(47r), where 7 is the Euler constant. 



Si 



6A 1 

r ~, B\ 



(4tt) 2 ?' 



6A 1 



4ff 2 1 
(4^) 2 e" 



4 3 2 1 



Ci 



(4tt) 2 e 
8 5 2 1 



(4tt) 2 e 



1 1 



647T 2 ? 



(4tt) 2 e 
,-D 2 ,i=3 



1 1 

64tt 2 ? 



£>3,1 = 



1 1 



64tt 2 



For all steps we always have 

A n = A 1 ,A n = A\,E n = Ei,B n = Bi,B n = Bi,F n = Fi,C n = Ci,Di : 2,3,n = D 1,2,3,1 
and for n > 1 

I„ = -1 + A U % = -l + B 1: W n = -l-Ci 



(107) 

(108) 
(109) 
(110) 

(111) 
(112) 



B. Comments related with the presence of the fermions in the game 



We must remark that, on eqs. (B4), (BE) and (B1C) of appendix H, the following terms 



&g 2 



,2 1 



8g 2 



■3m 



2 1 



"(47T) 



and 



2fcn 



and remarked by the authors 



i should be neglected. As stated by the authors of j2_ 

of ]12| , these terms will be canceled by contributions from higher order loops. Since we are concerned only about 
the one-loop approximation, we do not have to worry about them. Nevertheless, in the 0(4) linear sigma model, 
i.e., when g — 0, none of the above terms will be present, and our model will be order by order renormalizable in 
the regions of validity of the MSCR. This occurs because our tree-level resummed masses are related by a symmetry 
relation that always guarantees the cancellation of the UV divergences. 
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VII. NUMERICAL ANALYSIS 



In this section, we present numerical solutions of the gap equations for the tree-level meson and fermion masses 
and the condensate derived in section IV C including all diagrams which belong to the one loop order. 

As an approximation, only for the sake of obtaining continuous curves, in the numerical evaluation we considered 
All = also in the intermediate temperature region. Rigorously speaking, the curves should only be trusted in the 
low and high temperature regions. Figure || shows the tree- level resummed meson masses, eqs. (|79| ) and (|80|), as 
functions of the temperature. We show in Figure ^ the tree-level fermion resummed mass, eq. (|8l|), as a function 
of temperature. The tree-level masses behavior exhibit the fact that the MSCR has solved the problem of tachyonic 
masses. In Figure || the chiral condensate v, eq.([32|), as a function of temperature is shown whereas in Figure ^| the 
condensate is ploted in the case M n = 0. 

Since at low temperatures the condensate dominates, the mesons masses suffers its influence in this region. The 
sigma mass decreases and they approach each other to become degenerate in a temperature of about SOOMeV. This 
confirms the results we found in a phenomenological approach to the linear sigma model p0| ] . 

The condensate is a slowly decreasing function of the temperature, which is a signature of the order parameter when 
the symmetry breaking term is present. The qualitative behavior of the results shown in Figs. ^|and ||can be compared 
with the ones obtained by Chiku and Hatsuda jljj since OPT also sums three-point vertex diagrams, as our method 
does. Some differences may be attributed to the incorporation of the fermions, as performed in our method. In the 
absence of the chiral symmetry breaking term, i.e., when c = 0, the non- vanishing solutions of the extremum condition, 
eq. (H), are obtained numerically by equation @ with M% = 0, = 2\v 2 and M$ = M^,(M% = 0,M% = 2\v 2 ). 
The solution is depicted in Fig. |9| and gives an indication of first order phase transition. This result agrees with the 
predictions of first order phase transition found in previous analysis by Roh and Matsui |l8| , Petropoulos |l!| , Chiku 
and Hatsuda pl| , Randrup |3l]] and Bilic J|^] . Of course we have to bear in mind that our result is at one loop order 
in the perturbative expansion. It may well be that near the critical temperature higher order corrections become 
crucial and may change the order of the phase transition. 

The tree fermion mass, Fig. |?], does not become zero when chiral symmetry is restored and v ~ * since we 
considered contributions from the mesons, given by eq. (|o"o|). On the contrary, when the temperature is > 200MeV 
these contributions dominate the variable v and the fermion mass increases with temperature. The behavior of the 
fermion mass is in agreement with the results found by Panda in p3| for the quark meson coupling model. 



VIII. CONCLUDING REMARKS 



In this paper, we presented a modified self-consistent resummation (MSCR) at finite temperature. Results for the 
chiral fermion meson model and the massless A</> 4 model in the weak coupling limit were obtained and analyzed. We 
have shown that our procedure properly resumes higher order terms which cures the problem of the breakdown of the 
perturbative expansion. 

We have also shown that the MSCR, when applied to the study of the chiral fermion meson model, has the essential 
features which leads to the satisfaction of Goldstone's theorem and renormalization of the UV divergences, in the 
low and high temperature regions. We have explicitly shown that the scheme breaks down around T c i.e., in the 
region of intermediate temperatures. The application of the MSCR in these three physically different regions (low, 
intermediate and high temperatures) revealed a source of mistakes usually found in the literature, that is to treat all 
ranges of temperatures in the same way. It is naive to expect that the same approximations which is valid, e.g., for 
high temperatures would be enough in the intermediate temperature region, since quantum fluctuations are known 
to play a major role there. 

This division was essential to identify the regions where higher order terms and resummation are crucial. It is valid 
to remember that even when higher order loops are taken into account, the resummation is still necessary since the 
tree-level masses will become tachyonic even below the critical temperature (in theories with spontaneous symmetry 
breaking) and break the perturbative expansion. This breakdown of the perturbative expansion can also happen in 
massless field theories, like QCD, due to the appearance of infrared divergences. As we discussed, the breakdown of 
perturbative expansion in finite temperature field theory requires resummation techniques as the MSCR to recover 
the reliability of perturbative expansion. 

In each region renormalization and satisfaction of Goldstone's theorem were discussed in detail. In our study, we 
have also addressed a usually avoided point: the inclusion of the fermions. Finally, we have re-examined the chiral 
phase transition in static equilibrium in terms of the linear sigma model with our MSCR. 
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The gap equations for the tree-level masses, are constructed by our method and in the effective Lagrangian they are 
renormalized. For the particular case of intermediate temperatures region, the gap equations would be renormalized 
in the (reorganized) effective Lagrangian only if All = Ii a — H n = 0. In most of the approximations found in the 
literature, the gap equations are reached by some technique or via some ad-hoc procedure but the Lagrangian is yet 
the original one. This makes the renormalization process non-trivial, unless a finite cut-off is used and the theory is 
treated as an effective model |27],^] . As pointed out by Chiku and Hatsuda |Q , the resummation must be done also 
in the counterterms, which is essential to show the renormalization. 

At this point, it is extremely worth emphasizing that, although one has the freedom of adding and subtracting mass 
parameters to the Lagrangian, in this case they can not be completely arbitrary. If the mass parameters introduced 
were different for the pion and sigma fields (i.e., \M\a 1 + \Mi-k 1 and, of course, the same quantities subtracted, with 
Mi j^z M2), neither the 0(4) linear a model is renormalizable in any given order nor Goldstones theorem is satisfied. 
This will happen even if the mass parameters are determined by some physical condition as FAC or principle of minimal 
sensitivity (PMS). So, the most important fact behind the fulfillment of Goldstone's theorem and renormalizability 
of theories with SSB is the chiral symmetry that must dictate which mass parameter should be introduced to the 
Lagrangian. 
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As mentioned earlier the vacuum contribution to Qi(T, v) is divergent and requires renormalization. In this sub- 
section we also use dimensional regularization in the computation of the effective potential. 
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APPENDIX A: RENORMALIZATION OF THE EFFECTIVE POTENTIAL 




(Al) 




(A2) 



where L(m) is the usual zero temperature loop integral 




(A3) 



with d 4 p = dp4(Pp being the four Euclidean momentum. 

The divergent integral L(m) can be evaluated in the standard manner 




(A4) 



The quantity £l\(m) is then obtained with the integration of mL (m) 




(A5) 



With this expression we can find the zero temperature effective potential, 



18 



-) + 3(m CT <-> TOtt) - 8(m a <-> rrty) 



(A6) 



The renormalization of the thermodynamical potential at the end amounts to the determination of the parameters 
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APPENDIX B: ONE-LOOP SELF-ENERGY AT FINITE TEMPERATURE 



At zero momentum the expressions for the self-energies are given by 
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with Z defined as 
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FIG. 1. Figure 1: The logarithm of the interaction partition function 1PL Solid line is the tree-level pion meson propagator, 
the dashed line is the tree-level sigma meson propagator and the solid line with a arrow is the tree-level fermion propagator. 



23 





FIG. 2. Figure 2: One-loop self-energy for the pions. 
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FIG. 3. Figure 3: One- loop self-energy for the meson sigma. 
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FIG. 4. Figure 4: One-loop self-energy for the fermions. Note that the lowest corrections for the fermions are of order 
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FIG. 5. Figure 5: (a) The 1PI one-loop self-energy diagram of the X(f> 4 model and (b) a "daisy" 
attached bubbles which contributes to the self-energy. 
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FIG. 6. Figure 6: Tree-level resumed meson masses. 
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FIG. 7. Figure 7: Tree-level resumed fermion mass. 
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FIG. 8. Figure 8: Condensate v as a function of the temperature. 
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FIG. 9. Figure 9: Condensate v as a function of the temperature in the chiral limit (M^ — 
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